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ABSTRACT 

Large solid propellant rocket motor, meant for next generation launch vehicle was static tested at SDSC 
SHAR. During burning of 103 sec, 207 tonnes of HTBP based solid propellant motor was burnt and 
major exhaust pollutants are Aluminum Oxide (A1 2 0 3 ), Hydrogen Chloride (HC1) and Carbon Monoxide 
(CO). Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model is used to predict the 
dispersion of exhaust pollutant in the form of vapour and ground level concentrations. HYSPLIT model is 
run through a Graphical User Interface (GUI) for ease in setting up trajectory, concentration or deposition 
guidance. HYSPLIT was configured to simulate the pollutant dispersion from Rocket static test bed 
source location using meteorological input from the Meso-scale Model version 5 (MM5). HYSPLIT 
simulated from 08:00 1ST to 10:00 1ST on the day of test at ground level, 10m, 100m, 300m and 500m 
altitudes. GUI converted the MM5 meteorological grids into HYSPLIT binary format, model was run and 
post-processed the data to create output graphics and displayed in KMZ Goggle Earth. MM5 model 
initialized with local observations like GPS RS, AWS and MET tower observations and are run at 3 km 
spatial resolution and with thirty-one vertical levels. MM5 provided input into HYSPLIT that may better 
predict the unique meso-scale environment. Rocket motor emission rate, composition and duration inputs 
are configured in HYSPLIT. Attempts were made to simulate HC1, CO concentrations at elevations of 
10m, 100m, 300m and 500m altitudes and A1 2 0 3 particle deposition at ground. In future work, HYSPLIT 
dispersion simulations will be validated with the observations. 

1.0 Introduction 

Air quality monitoring is a paramount requirement of safety, health and environment divisions of SDSC 
SHAR, Andhra Pradesh, India. Since 2000, SDSC SHAR has constantly analyzing the atmospheric 
dispersion models for predicting the behavior of rocket exhausts in the atmosphere and propellant vapor 
dispersion from storages and as well as launch vehicles. The pollutant dispersion analysis are routinely 
used to assess the environmental impact of exhaust pollutants from rockets with respect to permissible 
toxicity thresholds and air quality standards. Beyond analyzing the environmental impact, the dispersion 
data is used to map the safe corridor for human life and work facilities. This paper deals with the exhaust 
dispersion analysis from large solid propellant rocket motor static test firing on 24 th January 2010 02:30 
UTC or 08:00 1ST (Fig.l) using HYSPLIT model over SDSC SHAR. 

2.0 Large Solid Propellant Rocket Motor (S200) 

The large solid propellant rocket motor, S200 consists 
of HTPB based composite propellant with 86% solid 
ingredient loading. There are three segments in the 
motor. The overall length is 21.9m. Diameter and 
propellant mass is 3200 mm and 207 tones respectively. 

When this solid motor of S200 rocket is static test fired, 
the major exhaust pollutants are Aluminium Oxide 
(A1 2 0 3 ) - 34%, Carbon Monoxide (CO) - 26.6% and 
Hydrogen Chloride (HC1) - 20.7%. The emission rate 
of the exhaust pollutant is considered as 2009 kg/s. 

Based on the above data the vapour dispersion of HC1, 

CO and particle deposition of A1 2 0 3 was analysed using 
HYSPLIT model at different altitudes. 



Fig.l View of S200 static test and plume rising 







3.0 HYSPLIT Model Description 

The HYSPLIT model has been developed at the NOAA Air Resources Laboratory (ARL). The model is a 
complete system for computing atmospheric trajectories, complex dispersion and simulations using either 
puff or Lagrangian particles (Draxler 2011). HYSPLIT applications include tracking and forecasting the 
release of pollutant material, air-parcel trajectories, volcanic ash, and smoke from wildfires. HYSPLIT 
uses gridded meteorological data to drive the trajectories and concentrations either forwards or backwards 
in time at regular time intervals. The HYSPLIT model uses grid meteorological data and dispersion rate is 
calculated from the vertical diffusivity profile, wind shear and horizontal deformation of the wind field. 
Air concentrations are calculated at the specific grid points for puffs and cell average concentrations for 
particles. The software is written mainly in Fortran with several modules in C code. 


4.0 Meteorological Model (MM5) data for Hysplit 

The meteorological modelling system is based on the PSU/NCAR Mesoscale Model (MM5). The well 
known MM5 system is a parallelized, limited area, non-hydrostatic, terrain following and sigma- 
coordinate model designed to simulate or predict mesoscale atmospheric circulation (Dudhia, J. et al., 
2005). Pre-processors are used to define computing grids, to assimilate global scale hydrostatic models 
outputs (global analysis and forecasts), topography, land use data and thus prepare initialization and 
boundary conditions input files for MM5. For Hysplit, NCEP/GFS global data with 3 hours and 0.5° X 
0.5° resolutions are used. 


The simulations are performed in a nested domains mode, which allows to consider in a single run both 
large calculation domains and high resolutions. The inner domains at high resolution are defined into a 
larger domain at lower resolution (coarse domain). Simulations are also performed in a two-way mode: 
calculations for all domains are carried out simultaneously, so that boundary conditions are provided at 
each time step to inner domains at high resolution. This method also allows feedback on the coarse 
domain. The horizontal resolution for the outer domain is 9 km and Inner domain is over Southern 
Peninsular India covering SDSC SHAR region with 3 km resolution. Both domains have 31 vertical 
levels defined from soil to lOOOOPa (approximately 15000m). 

The MM5 parameterization used is as follows: Reisner Graupel (Reisner2) explicit moisture scheme, 
Kain-Fritsch-2 cumulus parameterization, MRF planetary boundary layer scheme, RRTM radiation 
scheme, and five-layer soil model. MM5 model integrated on 24th Jan 2010 with 00UTC initial 
conditions, assimilation with GPS RS, meso- network AWS and other available observations. Model 
output hourly dumping made available to HYSPLIT runs for Static Test time at 02:30 UTC. Continues 
Meteorological Surface and soundings are measured from 00UTC on the test day. MM5 model derived 
Temperature, Relative humidity and wind profiles are validated with GPS-Radiosonde profiles shown in 
Fig.2 & 3. Validation results are confirming that MM5 model able to predict the lower tropospheric 
moisture intrusion and lower level strong easterlies due to Easterly wave activity over SHAR region. 


Comparison of Vertical profile of temperature 24-01-10 at 0130 UTC with GPS 



Comparison of Vertical profile of RH 24-01-10 at 0130 UTC with GPS 



Fig.2 Comparison of MM5 model derived & GPS-Radiosonde Temperature and RH profiles over SDSC 
SHAR region at 01:30 UTC 
























































Comparison of Vertical profile of Wind Speed 24-01-10 at 0130 UTC with GPS 


Comparison of Vertical profile of Wind Direction 24-01-10 at 0130 UTC with GPS 




Wind Speed (m/s) Wind Direction (Deg 

Fig.3 GPS RS Wind speed and direction profiles comparison with MM5 derived wind profiles for SDSC 

SHAR region at 01:30 UTC 


5.0 HYSPLIT Configuration 

Around 207 tones of solid propellant S200 rocket motor is statically test fired at ground level in 103 
seconds. MM5 model meteorological gridded data to HYSPLIT, ARL conversion program used to 
convert model required horizontal grid format. Configuration of HYSPLIT model is based on actual data 
and shown given below, the corresponding configuration of source emission rate, location, duration, setup 
and turbulence in HYSPLIT software is shown in Fig 4 and Fig 5. . 

Total Emission Rate : 60,00,000 kg/hr 
Total Heat Release : 10973765.8 W 

Release Duration : ~2 Min 

Pollutants: Rocket exhaust propellants are configured based on lab results, shown below 


Particles 

A1 2 0 3 (34%) 

Emission rate 

Diameter 

Density 

Molecular weight 

20,40000 kg/hr 
100 (am 

4 g/cc 

101.9 g/mole 

Gases 

HC1 (20.7%) 

Emission rate 

Density 

Molecular weight 

12,42000 kg/hr 
1.27 g/cc 

36.46 g/mole 

CO (26.64%) 

Emission rate 

Density 

Molecular weight 

15,98400 kg/hr 
1.84 g/cc 

28.01 g/mole 



Fig A HYSPLIT Configuration - source emission rate, location and duration 
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Fig.5 Setup & Turbulence - HYSPLIT Configurations 

6.0 Synoptic Weather on the day of static test 

On S200 rocket motor static test day, Easterly wave active over Southern peninsular India and effecting 
Coastal Tamilnadu and extreme South Coastal Andhra. In Fig.6, Kalpana-1 EIR image shows low clouds 
and 500m level Easterly wave flow over SDSC SHAR region. Static test time around 2:30 UTC, cloudy 
sky and isolated drizzles experienced over the region. 



Fig.6 Kalpana-1 EIR Brightness Temperature at 02:30 UTC and 24th Jan 2010 and MM5 Model analysis 

- 500m altitude mean wind flow 












































































7.0 Results and Discussion 

The firing of the solid rocket motor during the static test gives a large cloud of hot and buoyant exhaust 
products near ground level initially (Marco Cencetti, et al., 2007). The cloud rises and entrains ambient 
air until the cloud reaches an approximate equilibrium with ambient conditions like temperature, density, 
speed and turbulence. The cloud is then called the stabilized cloud. Based on Morten, Taylor and Turner 
model, stabilized cloud base altitude or the plume rise is calculated as given below: 

Plume Rise H = 2.675 Q 1/4 

Q = m Cp (Te - Ta) 

Where 

H = Final height in metre 
Q = Rate of release of heat in joules 
m = Emission rate of exhaust gas in kg/s 
C p = Specific heat of exhaust gas in J/kgK 
T e = Exhaust gas temperature in K 
T a = Ambient temperature in K 

Based on the above model, plume rise achieved was around 1238.4 m and assumed that stabilized cloud 
base is around 600 m altitude. The HYSPLIT simulation results depicted in Fig 6,7 and 8. 



Fig.6 Deposition of Al 2 Q 3 at ground level, 10m, 300 m and 500m altitudes 




















HYSPLIT simulation of A1 2 0 3 pollutant deposition in 5 km range of test bed is about 1.2 g/m 2 . From test 
bed in 20 km range deposition is about 1.0 g/m 2 . A1 2 0 3 concentrations at 10m altitude is about 0.1 mg/m 3 
up to 10 km range, at 300m altitude is about 0.017 mg/m 3 up to 10 km range and at 500m altitude is about 
0.01 mg/m 3 up to 10 km range. The simulation results of deposition and dispersion of AL 2 0 3 pollutant is 
uniformly distributed from surface to 500 m altitude. 

HC1 concentration at 10 m altitude simulation explains southward movement of plume with maximum 
concentration of about 1.4 ppm within 5 km range. At 500m altitude HC1 concentration simulation shows 
the westward movement of the exhaust plume. The dispersion of the exhaust plume at upper level (500m 
altitude) is following the mean flow of the easterly wind, which was confirmed in local wind observation 
and MM5 simulations. 




Fig. 7 Dispersion concentration ofHCl at 10m and 500m altitudes 


10m and 500m altitude concentrations of CO within 5 km are 2.35ppm and 0.25ppm respectively. In this 
case also upper level concentrations is followed the mean wind flow direction is shown in Fig 7. 




Fig. 8 Dispersion concentration of CO at 10m and 500m altitudes 


8.0 Conclusion and Future Works 

The configuration of the HYSPLIT model is very generic, it could be set up to perform a variety 
of different scenarios. In general Lagrangian models are well suited for quick calculations from 
pollutant point sources. The dispersion of A1 2 0 3 , HC1 and CO of S200 solid propellant rocket 
motor firing at different altitudes from surface to 500m is simulated and results are presented. 
The comprehensive validation of the HYSPLIT would be taken up with experimental 
observations for the further understanding of the short-range dispersion and model validation. 
























This may pave the way to stride towards the prediction of exhaust gas behavior of largest solid 
motor and future launch vehicles of Indian Space Research Organisation. 
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